Quasiclassical approach to the spin-Hall effect in the two-dimensional electron gas 
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We study the spin-charge coupled transport in a two-dimensional electron system using the 
method of quasiclassical (^-integrated) Green's functions. In particular we derive the Eilenberger 
equation in the presence of a generic spin-orbit field. The method allows us to study spin and 
charge transport from ballistic to diffusive regimes and continuity equations for spin and charge are 
automatically incorporated. In the clean limit we establish the connection between the spin-Hall 
conductivity and the Berry phase in momentum space. For finite systems we solve the Eilenberger 
equation numerically for the special case of the Rashba spin-orbit coupling and a two-terminal ge- 
ometry. In particular, we calculate explicitly the spin-Hall induced spin polarization in the corners, 
predicted by Mishchenko et al. Il3l . Furthermore we find universal spin currents in the short-time 
dynamics after switching on the voltage across the sample, and calculate the corresponding spin-Hall 
polarization at the edges. Where available, we find perfect agreement with analytical results. 



I. INTRODUCTION 



In the presence of spin-orbit coupling spin currents and 
spin polarizatio n ca n be generated as a response to elec- 
tric fields 1.2.3 4 5]. Recently the spin-Hall effect, i.e. a 
spin current that flows perpendicular to an applied elec- 
tric field has been observed experimentally in electron 
doped semiconductors and in a two-dimensional hole sys- 
tem . Theoretically one may distinguish the extrinsic 
from the intrinsic spin-Hall effect, depending on whether 
spin-orbit coupling arises due to scattering by impurities 
or from the intrinsic band structure of the samples. The 
intrinsic spin-Hall effect has first been studied for holes in 
p-type semiconductors in £| and for electrons in n-type 
semiconductors in [jjj]. In both cases, the striking result 
is the independence of the spin-Hall conductivity o s h of 
the strength of the spin-orbit coupling, at least when dis- 
order effects are ignored. Sinova et al. @ found in the 
two-dimensional electron system with Rashba spin-orbit 
coupling a universal value for the spin-Hall conductiv- 
ity, u s h — e/87r. Soon it was realized that the univer- 
sal spin-Hall conductivity still exists in the presence of a 
Dresselhaus term , and even in the presence of a weak 
in-plane magnetic field [ljj ■ ^ was pointed out that this 
may be related to a Berry phase in momentum space, i.e. 
the winding number of the spin-orbit field when going 
once around the Fermi surface. Shytov et al. showed 
such a connection explicitly in the specific case where the 
modulus of the spin-orbit field remains constant on the 
Fermi surface. 

Clearly it is an important question to ask how u s h 
depends on disorder. For the Rashba model the effect 
is quite dramatic, namely an arbitrarily weak amount 
of disorder fully sup presses the spin-Hall conductivity 
|!2ll3ll4ll5llrJl7llil| . Meanwhile it is understood that 
this surprising result is a special property of the Rashba 
Hamiltonian and is related to the linear-in-momentum 
spin-orbit field: The time derivative of the total spin is 



proportional to the spin c urrent , so that in a steady state 
both quantities are zero [HjlSl- In the case of a more 
general spin-orbit field, a finite spin-Hall effec t has been 
reported even in the presence of disorder [Illl8l2fi21l22) . 

With the conventional definition of the spin current, 
given by the anticommutator of the velocity operator 
and the Pauli matrices, the spin current is not conserved. 
Hence it is not automatically guaranteed that a current 
in the bulk induces a polarization at the edges of the 
sample. Furthermore the spin current is not directly ac- 
cessible experimentally; the measurable quantity is the 
spin polarization instead. Therefore it is of interest to 
study directly the electric field induced spin-density. A 
strategy followed by some authors is to discretize the 
Rashba Hamiltonian in terms of a tight-binding model 
which is then studied near the band edge. In this way 
the spin-Hall induced spin accumulation in systems with 
linear dimensions of several tens of the Fermi wavelength 
have been studied |23u24l ] . Macroscopic systems are con- 
venie ntly described in te rms of semiclassical kinetic equa- 
tions |llll3ll4l2ll22l25| . and in particular in terms of 
diffusion equations describing the coupled dynamics of 
spin and charge degrees of freedom. In clean systems, 
however, the spin relaxation length can become compa- 
rable to the elastic mean free path. In this situation the 
spin dynamics is not diffusive, and one has to go beyond 
the diffusive approximation. In what follows we present 
a theory of the spin-Hall effect in terms of quasiclassi- 
cal Green's functions, which covers the full range from 
clean systems to the diffusive limit. The first step is 
the derivation of the equation-of-motion, the Eilenberger 
equation. From there we derive the continuity equation 
for charge and spin, and obtain explicit expressions for 
the current densities. For the Rashba model it follows 
that no spin current can flow in a time and space inde- 
pendent situation. In the clean limit we find a solution 
of the equation-of-motion which corresponds to the uni- 
versal spin-Hall effect. In the general case of disordered 
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finite systems we solve the equation-of-motion numeri- 
cally. 

II. THE EILENBERGER EQUATION 



We evaluate the ^-integral explicitly in the limit where 
|b| is small compared to the Fermi energy. Since the 
main contributions to the ^-integral are from the region 
near zero, it is justified to expand b for small £, b w 
bo + £<%bo, with the final result 



We start from the Hamiltonian 

„2 



H = 



P 
2m 



(1) 



where b is the internal magnetic field due to the spin- 
orbit coupling and u is the vector of Pauli matrices. In 
the Rashba model for example b = ap x e z . For a spin 
1/2 particle one can write the spectral decomposition of 
the Hamiltonian in the form 



H 



-)(- 



(2) 



where e± = p 2 /2m± |b| are the eigenenergies correspond- 
ing to the projectors 



l±)(±l 



1 ±b ■ er 



(3) 



where b is the unit vector in the b direction. We write 
Green's functions in Wigner coordinates, G = G(p,x), 
where p is the Fourier transform of the relative coor- 
dinate and x is the center-of-mass coordinate. For the 
Green functions we make the ansatz 



G 



G R 





QK 



l_r 





— u 



g A 



where the curly brackets denote the anticommutator. 
Gq' A are retarded and advanced Green's functions in the 
absence of external perturbations, 



G 



R(A) 



1 



e + n - p 2 /2m b a Z R ( A ) ' 



(5) 



and T, R ^ are the retarded and advanced self-energies 
which will be specified below. The ansatz guarantees 
that in equilibrium the matrix of Green's functions with 
small letters is 



9 = 



1 2tanh(e/2T) 
-1 



(6) 



The main assumption for the following is that we can 
determine g such that it does not depend on the modulus 
of the momentum p but only on the direction p. Under 
this condition g is directly related to the ^-integrated 
Green function which we denote by g, 



9 = ~ d£G, i = V A l2m~ [ i. 

7T 



(7) 



For convenience we suppressed in the equations above 
spin and time arguments of the Green function, g — 
9tis 1 ,t 2 s 2 (P! x )- I n some cases Wigner coordinates for the 
time arguments are more convenient, g — > g SlS2 (p, e; x, t). 
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{1 - <9 e b • cr,g} 
{l + <%b -<r,<7}. 



(8) 
(9) 



In the equation of motion we will also have to evaluate 
integrals of a function of p and a Green's function. As- 
suming again that |b| <ti tp we find 



- / d£/(p)G«/(p + )ff + + /(p_)g_ 



(10) 



where p± is the Fermi momentum in the ±-subband 
including corrections due to the internal field, |p±| ~ 
p F T \b\/v F , and 



■ 9+ 



(11) 



Following the conventional procedure^ we derive now 
the equation-of-motion for g. From the Dyson equation 
and after a gradient expansion the equation-of-motion for 
the Green function G reads 



d t G+-{ - + —(b-(T),—G 
2 (m op ax 

= -ip,G]. 



i [b ■ cr, G] 



(12) 



The Boltzmann equation or Boltzmann-like kinetic equa- 
tions are obtained by either integrating (|12|) over energy 
e or over £, see |2^. Ref. for instance follows the first 
route, whereas we integrate over £. Retaining terms up 
to first order in |b|/ei? leads to an Eilenberger equation 
of the form 



o \ — ^ Try " ' a h -zr9v 

2 ^ to ap ax 



i[b w • cr,^]) 



(13) 



In the entire article we will take the self-energy as S = 
— i(g) /2t , which corresponds to s-wave impurity scatter- 
ing in the Born approximation; (...) denotes the angular 
average over p. 

To check the consistency of the equation we study at 
first its retarded component in order to ver ify that g R = 1 
solves the generalized Eilenberger equation. From Eq. (JSJ) 
we find that g R ' = 1 — 9^(b ■ cr) and using (fTTf) we arrive 
at 



g« = (lTdsb)[-±h ± -a 



(14) 
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Apparently both the commutators on the left and on the 
right hand side of the Eilenberger equation are zero, at 
least to first order in the small parameter bo. Similar 
arguments may also be used to verify that the equilib- 
rium Keldysh component of the Green function, g K — 
ta,nh(e/2T)(g R — g A ), solves the equation of motion. 

In the appendices we demonstrate how the frequency 
dependent spin-Hall conductivity and the equation-of- 
motion in the diffusive limit can be obtained from 
Eq. l|13l) . For the Rashba model our results agree with 
Ref. [13 . 



III. CONTINUITY EQUATION - VANISHING 
SPIN-HALL CURRENT 

Equipped with the Eilenberger equation it is not dif- 
ficult to see that for a spin-orbit field of the Rashba or 
linear Dresselhaus model the spin-Hall conductivity is 
zero. The argument is analogous to that of Ref. and 
makes use of the continuity equation. 

When taking the angular average of the Eilenberger 
equation (|13|) . the term on the right hand side vanishes 
and we are left with a set of continuity equations for the 
charge and spin components of the Green function. With 
9ss' = galas' + g ' fss' the equations read 

d t (go)+d x -J c = (15) 
d t {g x ) + d x -J* = 2^(b„xg„), (16) 

v=± 



dt(g-y) +d x -jy = 2^(b„xg„) 9 (17) 

v=± 

d t (g z ) + d x ■ J z s = 2^(b„xg„) z (18) 



with 



E 

v=± 



1 \Vv 



d 



2 \ m dp 



(19) 



The densities and currents are related to the Keldysh 
components of (g) and of J c jS integrated over e. Explic- 
itly the particle and spin current densities are given by 

— Jf(e;x,t) (20) 

= -~^ J ^Jf( £ ;x,t) (21) 

with No = m/27r being the density of states of the two- 
dimensional electron gas. In the the absence of spin-orbit 
coupling (b = 0) one recovers the well known expressions 



jo(x,t) = de(v F g«) 
jl(x,t) = —No J de(v F g?) 



(22) 
(23) 



In the presence of the field b the expressions are in gen- 
eral more complex. For the Rashba model, for example, 
the particle current is given by the lengthy expression 

jc(x,t) = ~N J de[v F {pgK) 

+a(e z x ( S K ) - (P(P ■ e z x g K ))}. (24) 

Finally let us consider the spin current with polariza- 
tion in z-direction for the model with both a Rashba and 
a Dresselhaus term for which the spin-orbit field reads 




Py \ I Px 

" | ~Px +0 -Py 

o V o 



(25) 



Because the field lies in the x— j/-plane, 3 Z S is simply given 
by Jf = (vfIJz)- Besides the field is just linear in p. As 
a result we find that the source term in the continuity 
equations (|16fl and l|17(l can be expressed in terms of the 
spin current: 

HE) = 2(b yfi g z ) = -2maJ z s x - 2mpJ z y (26) 
113 = -2(b xfi g z ) = -2ma.Jl y - 2m/3J* x . (27) 

In a stationary situation and for a spatially homogeneous 
system the left hand side of the continuity equation is 
zero and this implies a vanishing spin current, 



3a 



: 3 a 



= 0. 



(28) 



IV. CLEAN LIMIT - UNIVERSAL SPIN-HALL 
CURRENTS 

We consider now the Eilenberger equation in the clean 
limit, t — > oo, and study the linear response to an ho- 
mogeneous electric field. For a realistic system with at 
least weak disorder this study still gives reliable results 
on short time scales, i<r. Generally an electric field 
can be included in the quasiclassical equations-of-motion 
by the substitution <9 X — > d x — \e\Ed e . The Keldysh com- 
ponent of the linearized Eilenberger equation becomes 



v=± 

2 



Y,(dt9?--V-Pvd e g?™ (29) 
z — ' m 

{(E ■ d p )(b v ■ cr),d e g^} + i[K ■ <r,g«}) = 0. 



In the following we focus on the spin components of the 
equation. Explicitly we get 



9t9x = 2b v,og 



K 



HE- 



Pb xfi - v F dtb x ,o + d p b Xj o F t (30) 



dtg y = -Zbxpg 
+|e|E. 



A" 



P&y,o - v F dt:b y fi + dpb yt o F e (31) 



d t gf = 2(b x ^~b y ^) 



+2(b X fid(by t o - b v fld^b Xt o)go- 



(32) 



where for the sake of brevity P = ^ vp v /2m and F e 
2d e tanh(e/2T). For the g^ component one obtains 



2 n K 



2„K 



dt 2 



2F e \e\ [6 Xi0 (E • d p )b yfi - b yfi {V ■ d p )b x , ] ■ 

(33) 

Notice that only the second of the two terms involving 
the electric field in Ea. (|29f) remains in the equation for 
the g% component. The solution of this differential equa- 
tion is the sum of an oscillating and a time independent 
term. Due to the (undamped) oscillations it is clear that 
a stationary solution is never reached so the arguments of 
the previous section leading to vanishing spin-Hall cur- 
rent do not apply. The time independent solution of the 
differential equation is related to a zero-frequency spin 
current given by, 



js = -^(PF(E-a p )vi>), 



tan* = b yfi /b xfi . (34) 



Notice that the spin current does not depend on the mag- 
nitude of the field b, but only on the variation of its di- 
rection when going around the Fermi surface. An even 
more explicit result is obtained when the spin-Hall con- 
ductivity tensor is antisymmetric 



1 / y,x x,y\ 
°sH = -Kjj -<7 S S) 

\ e \ 

= --^((PFydp, -pFxd Py )^) 

8tt/ 2tt p ' 



(35) 
(36) 
(37) 



i.e. the spin-Hall conductivity is the universal number 
|e|/87r times the winding number of the internal field b 
when going once around the Fermi surface. 

We notice that Eq. is consistent with 8.9] and 
also with |lOj | where the spin-Hall conductivity ignoring 
disorder has been calculated using the Kubo formula for 
a Rashba-Dresselhaus system in the presence of an in- 
plane magnetic field. Eq. (|37l) which relates the spin-Hall 
conductivity with a winding number, i.e. the Berry phase 
in momentum space, generalizes the equivalent result of 
[TH ] , where it has been assumed that the modulus of b is 
constant on the Fermi surface. 



V. FINITE AND DISORDERED SYSTEMS - 
NUMERICAL RESULTS 

In this section we solve Eq. H13fl numerically for 
the Rashba model. Compared to the diffusion 
equationsi^*2iiS^, one advantage of our method is that 
we have access to length scales that are shorter than the 
mean free path I = vft. This is crucial in weakly disor- 
dered systems, |b| > 1/t, where the characteristic length 
scale for the spin polarization, the spin relaxation length, 
is of the order of the mean free path. Furthermore, when 




FIG. 1: The two-terminal geometry under consideration: A 
rectangular strip of a two-dimensional electron gas is con- 
nected to two reservoirs. We assume that the strip and the 
reservoirs are made of the same material, i.e. the spin-orbit 
field exists also in the reservoirs. 



considering time-dependent situations we can study the 
time evolution on time scales which are shorter than the 
scattering time r. 

In the following we will consider a geometry as shown 
in Fig.^ A rectangular strip of length L x and width L yi 
is connected to leads at x = and x — L. At interfaces 
the Eilenberger equation has to be complemented with 
boundary conditions22i22i2S. The boundary condition be- 
tween the strip and the leads is obtained assuming that 
the leads are made of the same material as the strip, i.e. 
there is no Fermi surface mismatch, and that both leads 
are in thermal equilibrium. For directions p pointing into 
the strip the Green function at the interface (x = 0, L) 
reads 



g K (p in :x.) 



=0,L 



5™(pm;x) 



x in the lead 



(38) 



where V is the applied voltage. At the boundary with 
an insulator the Green functions for in and outgoing di- 
rections are related via a surface scattering matrix, 



Sg K ( Pin )S' 



(40) 



Generally, the S'-matrix can be calculated by solving 
the quantum mechanical surface scattering problem. For 
simplicity we assume specular scattering and assume that 
boundary scattering does not induce transitions between 
the two spin-orbit subbands, so that 



Ik in ±) 



exp(±h9)|k out ±) 



(41) 



as it is found for smooth confining potentials J30j. Ex- 
plicitly the surface S'-matrix for the Rashba Hamiltonian 
is 



S 



e lv cos $ — sin i? 
sin^ e^^cosf? 



(42) 



where the angle ip characterizes the ingoing direction, 
Pin = (cost/?, sin ip). Our numerical results are obtained 
assuming that the relative phase-shift in the scattering 
for the two subbands is negligible, i.e. ?9 = 0. 
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FIG. 2: Spin polarization in the presence of an electrical cur- 
rent flowing in ^-direction for a strip of length L x = 20/ and 
L y = 10/. The spin-orbit coupling strength is a = 10 _3 «f 
and the elastic scattering rate is 1/r = qpf/2. The spin po- 
larization is given in units of the bulk value, So = —\e\EaTN . 



Finally, to integrate the equation of motion numeri- 
cally we have to discretize the space coordinate x and the 
Fermi surface. In dirty systems g K (p) is nearly isotropic, 
so it is clear that a few discrete points pi on the Fermi 
surface are sufficient. In clean systems this is not a priori 
evident, but numerical tests show that even in this case 
convergence is reached quickly. Typically we describe the 
Fermi surface with a set of twenty to forty p^. 

First we show numerical results for the spin polariza- 



FIG. 3: S y in units of So as a function of x for L x = 200/, L y = 
100/, a/v F = 10~ 3 and ap F r = 0.005,0.01,0.02,0.05,0.1,1 
(from bottom to top). 



tion in the stationary limit. Fig. [3 depicts the voltage in- 
duced spin polarization for L x = 201, L y = 10/, cxpft = 2 
and a/vF = 10 -3 ; all our results are linear in the applied 
voltage, due to the linearity of the underlying equations. 
In the bulk, only the 5,, component is nonzero, and given 
by 5*0 = -\e\EarNo (afl^ . A spin-Hall effect induced 
spin polarization is found in the corners, as it is expected 
in 13] . The spin polarization however is not purely in 
z-direction but has also components in ir-direction. 

Fig. shows S y (x,y — Ly/2) by varying disorder. In 
the diffusive limit and assuming that the spin polariza- 
tion vanishes at the interface to the leads, it has been 
predicted that 0] 



S y (x) = S I 1 - 



cosh[(x - L x /2)/L s 
cosh(L x /2L s ) 



(43) 



where L s is the spin relaxation length. Apparently with 
the boundary condition we choose a spin polarization still 
exists near x — 0, L Xl in particular in the clean limit. 
Some mean free paths away from the interface on the 
other hand the data can be well fitted with an exponential 
increase or decrease, both in the clean and dirty limit. As 
a result we obtain the spin relaxation length as a function 
of disorder, shown in Fig.^J In the dirty limit, apFT <C 1, 
our numerical result agrees with what is expected from 
the diffusion equation, L s = \JDt s = l/2apFT. In the 
clean limit, for which we are not aware of any quantitative 
predictions, the spin relaxation length is of the order of 
the mean free path, L s w 1.27/. 

In the following we will consider the time evolution of 
the spin polarization and the spin current. We start with 
a system in thermal equilibrium, switch on the voltage 
and observe the relaxation of the system into its sta- 
tionary non-equilibrium state. It is a nontrivial problem 
to describe such a situation theoretically. One might be 
tempted to allow a time dependent voltage in the bound- 
ary condition, Eq. (|39|l , and to follow then the time evo- 
lution. In this case the charge density becomes time 
dependent and inhomogeneous. This procedure makes 
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FIG. 4: Spin relaxation length L s in units of I as a function 
of disorder, obtained by fitting the spatial dependence of the 
electric field induced spin polarization (shown in Fig. [3J| using 
S y = a + &[exp( — x/ L a ) + exp( — (L x — x)/L 3 )]. The diffusive 
limit expression is shown as a dashed line. 



sense for non-interacting electrons, but not for interact- 
ing electrons where the long range Coulomb interaction 
enforces charge neutrality. In principle, the interaction 
can be included into the quasiclassical formalism explic- 
itly, see e.gi^i, however this is beyond the scope of the 
present paper. 

Instead, we assume in the following that a voltage 
difference across the leads instantly results in a homo- 
geneous electric field in the sample. Thus one has to 
solve Eq. with the initial condition <?(p,e;x, t) — 
t&nh(e/2T)(g R — g A ) and taking into account the elec- 
tric field via the substitution <9 X — > d x — \e\Ed e . In the 
numerics, however we find it more convenient to work in 
a scalar gauge, since then the (static) electric field dis- 
appears from the equation-of-motion and is present only 
in the initial condition and in the boundary condition. 
Generally the gauge transformation for the fields and the 
Green function reads 



3*1*2 



A + 9 xX 

<t> - dtX 
exp{-i|e|[x(*i) 



(44) 
(45) 

X(t2)}}g tl t 2 - (46) 



In the end we have to solve Eq. with the boundary 
condition (|39() and the initial condition 



9 (P,e;x,t 



0)=tanh(i±MM) {9 «- 9 \ (47) 



where (f>(%) interpolates linearly between the two leads, 
<j>(x)=V(L x /2-x)/L x . 

In Fig. JSJ we show the spin current j' y as a function 
of time in the bulk and at the interface to the leads of 
a rather clean system (olpft — 2, a/vp — 10~ 3 ). On 
short time scales the bulk current agrees with what we 
found ignoring disorder in Sect. II VI The spin current os- 
cillates as a function of time with frequency 2apF , the 
time average is given by the universal spin-Hall conduc- 
tivity. In the bulk, for the weakly disordered system we 
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FIG. 5: Time evolution of the spin-Hall current at the in- 
terface to the leads and in the bulk. In the bulk we compare 
our numerical result (data points) with the analytical result 
(full line) of Eq. (1481 . Near the leads, only numerical data are 
available (dashed curve). j% <y is evaluated at y = L y /2, x — 
(boundary) and x — L x /2 (bulk) for L x — 201, L y — 101, 
a/vF = 1CF 3 and appr = 2. 



are considering, the time-dependent spin-current is given 
by 



3 



s.y 



\ejE 
8tt 



[exp(-t/2r) - exp(-3£/4r) cos(2ap F t)} 



which can be obtained from the frequency dependent 
spin-Hall conductivity given in the appendix. On the 
time scale of the spin relaxation time, here given by the 
scattering time r, the bulk spin current becomes expo- 
nentially suppressed and goes to zero in the stationary 
limit. Near the leads, on the other hand, the situation is 
somewhat different, since a finite spin current remains in 
the stationary limit. 

An important question is whether the spin current po- 
larizes the electron system at the edges. In Fig. [S] we 
show the spin polarization in z-direction across the sys- 
tem at x — L x /2 as a function of time. Since in the early 
time evolution spin current flows in the bulk, spin density 
accumulates near the edges. When the spin current dis- 
appears also the polarization vanishes. We see that the 
spin polarization at the edges oscillates as expected with 
frequency 2app- In the cleaner systems oscillations are 
of course faster. Remarkably, the maximum amplitude of 
oscillation relative to the bulk value is larger in the dirty 
system (apFT = 0.25 ), where it is almost of the order of 
one. 

This can be understood as follows: a rough estimate of 
the spin polarization at the edge is s z ~ T s jl y /L s . With 



t s - T/(ap F T) 2 
fs, v ~ eE(ap F T) 2 
L s - l/{ap F T) 



(49) 
(50) 
(51) 



the result is indeed s z ~ Sq = eEarNo. In the clean 
limit, on the other hand, the typical time and length 
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FIG. 6: Spin-Hall effect induced spin polarization S z in units 
of So as a function of y and t at x = L x /2 for L x = 20/, 
Z/y = 10/, cx/vf = 10~ 3 , and apFT = 0.25,2,5 (from bottom 
to top). 

scales are t s ~ r and L s ~ Z, from which we estimate s z ~ 
Sq/ {apFT), in agreement with our numerical findings. 

We expect the detailed structure of the spin polariza- 
tion near the edges to depend on the boundary condi- 
tions, i.e., on the type of the confinement at the edges. 
How boundary conditions affect the spin polarization at 
the edges is certainly an experimentally relevant issue, 
which however is beyond the scope of this paper. 

VI. SUMMARY 

We studied the spin-Hall effect in a two-dimensional 
electron system by applying the method of quasiclassical 
Green's functions. The method has its strength in the de- 
scription of macroscopic systems, i.e. systems with linear 
dimensions that are large compared to the Fermi wave- 
length. In particular we derived an Eilcnbcrger-like equa- 
tion in the presence of a generic spin-orbit coupling. We 
also showed that the method allows to derive in an ele- 



gant way various results present in the literature. For the 
special case of the Rashba model we calculated numer- 
ically the spin-Hall current and the spin-Hall effect in- 
duced spin polarization on a strip. From our data we were 
able to extract quantitatively the spin relaxation length 
in the entire regime from the clean to the dirty limit, 
which is not covered by the diffusion equation approach. 
Although in the Rashba model the zero-frequency spin- 
Hall conductivity is zero we found a spin-Hall induced 
spin polarization at the edges on a short time scale after 
switching on the voltage. We believe that our approach 
will stimulate further work in the field. 

We thank U. Eckern for valuable discussions. This 
work was supported by the Deutsche Forschungsgcmcin- 
schaft through SFB 484. 

APPENDIX A: DIFFUSIVE LIMIT 

When spatial and temporal variations are slow it is 
in many cases convenient to study the equations of mo- 
tion in the diffusive limit. Usually, in this limit the full 
angular dependent Green function g(p) can be easily con- 
structed from its angular average (<?(p)) so it is sufficient 
to study (<?(p)). To determine the equation for (<?(p)) 
we first write all the terms of Eq. with g on the left 
hand side and those depending on (g) on the right hand 
side. For the Keldysh component we hence get 

(M + Mi )g K = (N + Ni ) (g K ) ( Al ) 

where 

M g K = g K + T d t g K + v F Tp ■ 8 x g K 

+ ir[h -a,g K ] (A2) 

= -irj^lP-apCbo^),^} 
-hr [d $ (b - *), {b -cr,g K }] 
-i{(d 5 b .«x), 5 A '} (A3) 
^o<3 A ') = (g K ) (A4) 
Ni{g K ) = —tebo-o-, (A5) 

Here Mi and Ni are small in the expansion parameter 
|b|/ejr. The Eilenberger equation is then rewritten as 

g K = (Mo + Mi ) " 1 (N + Ni ) (g K ) , ( A6) 

i.e. to first order in |b|/ep 

g K = (AV 1 + M^N, - M - 1 M 1 M Q - 1 ) (g K ), (A7) 

from which the equation for the s-wave component of the 
Green function becomes 

(1 - (AV 1 ) - (M^Nx) + (M-^hM- 1 )) (g K ) = 0. 

(A8) 
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In the low frequency, long wavelength limit this is the 
generalized diffusion equation obtained in the literature 
in various limit giiSLSi. The explicit form is obtained by 
evaluating the angular average of the operator product 
M~ X N. 



with 



In the Rashba model for instance where b 
one finds 



Mo 



L 














L 





ap x 








L 


ap y 


- 


ap x 


-apy 


L 





Qy 


Qx 




Q u 








'•) 


Qx 



















0/ 



= ap x e 2 



(A9) 



L = 1 + rd t + vpTp ■ <9 X 
a = lappT 

Qx,y = aT(d x ,y - (p • d x )p x ,y) 



and 



Nn 



1 

-apy/v F 
apx/v F 




(All) 
(A12) 
(A13) 



-ap y /vF aftx/vF 

1 

1 

1 



(A10) In a dirty system, where a <C 1, the result reads 



(A14) 



dt - Dd 2 

-2Bdy d t 

2Bd x 




~2Bd y 
Ddl + 


2Cd x 



2Bd x 


2ca, 





-2Cd x 

-2ca 



y 



d t - Ddi + 2t: 



( G?o K > \ 

V (gh J 



= o 



(A15) 



where D = \v f t is the diffusion constant and 
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1 + a z 



v F a 



1 1 a 1 



2(1 + a 2 ) 2 ' t s 2rl + a 2 ' 



(A16) 



In the clean limit (a ~ 1) there is no spin diffusion, and 
Eq. (|A1 5|> is not justified. However the equation is con- 
structed in such a way that it applies in the clean limit 
for a space and time independent spin polarization. 

In the general case (arbitrary a) we investigate spin 
relaxation for a spatially homogeneous system, where the 
equation of motion reads 



L{1? 



a 2 )-L 2 --a 2 
' 2 



(A17) 



[L(L 2 +a 2 )-L 2 ](g?) = (A18) 



with L = 1 + rdt ■ Clearly the spin-dynamics for the each 
component is determined from three relaxation times. 
For the ^-component, e.g., these are 



1 

Tl 
1 



_L±_L^i_4 a 2. 

T 2 , 3 2r 2r 



(A19) 
(A20) 



For dirty systems the longest of these times is T3 w r/a 2 , 
in agreement with what we find from the diffusion equa- 
tion ljA15|l . In the clean limit the "relaxation" time be- 
comes complex corresponding to an oscillating spin po- 
larization. 



We go back to (|A1(I and solve it for an infinite system 
under the influence of a uniform but time-dependent elec- 
tric field in x-direction. We choose the vector gauge as 
in section HV1 i.e. <9 X — > <9 X — |e|J5e? £ . To linear order 
in the external field and transforming to Fourier space, 
(|A1|) becomes: 



M g K = (l + N 1 ){g K )+S E 



(Bl) 



where Mo and Ni are the same as before, Se is a source 
term due to the electric field: 



E 




(B2) 



with E = \e\lE(uj)d e gQ' cq . Inverting M and performing 
the angular average one obtains 



(Pv9 



K\ 



2{L 2 + a 2 ) 



v F 



(B3) 



where the two terms in square brackets correspond to 
bubble and vertex corrections, respectively, in the dia- 
grammatic language. Furthermore we see that the spin 
polarization along e y contributes to the spin-Hall current. 
The spin polarization along e y is obtained by inverting 
Mq and performing the angular average 



(Sf ) = 



a 2 E 



«f2[L(L 2 + a 2 )- (L 2 + f )] 



(B4) 
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In the static limit (L = 1) this reads (c/y) = ^E, 
corresponding to the bulk value So, quoted in the text. 
This value for (g*f) leads to the cancellation discovered 
previously^ and to the vanishing of the static spin Hall 
conductivity. Finally, by combining Eqs. 1B3IB4I) . one 
obtains the frequency-dependent spin-Hall current with 
polarization along e z as 



fs, y {u) = a sH (uj)E(uj) 



8n L(Li + a 2 ) - (L 2 + 



E{uj) (B5) 



which agrees with the result given by both Mishchcnko et 
alpLi and Chalaev and Lossit To examine the transient 
response in time, we use E(t) — Edit) so that 



(B6) 



where 



f s , y (t)=E dt'a sH (t') 



(B7) 



The poles of the integrand may be found by express- 
ing lot = x + iy and setting to zero the real and imagi- 
nary part of the cubic polynomial in the denominator of 

x[x 2 - 3y 2 - 4y - (1 + a 2 )} = (B8) 
x 2 (3y + 2) - y 3 - 2y 2 - y(l + a 2 ) - y = 0. (B9) 

In the large a limit, one obtains the three solutions 

x = 0, y = -1/2; x = ±a, y = -3/4. (BIO) 
Computing the residues with the same accuracy yields 



= _\e[ejti /2r 

8tt 2t L 
-2ae- 3t / 4r sin(at/T)]. 



<TsH(t) = -p^[e-*/2r_ e -3t/4r cos(at/r) 



(Bll) 



Finally, inserting the above result into Eq. (|B6|> . one ob- 
tains - to leading order in 1/a — the time-dependent cur- 
rent quoted in the text. 



1 M. I. Dyakonov and V. I. Perel, Phys. Lett. 35A, 459 
(1971). 

2 V. M. Edelstein, Sol. Stat. Commun. 73, 233 (1990); J. 
Phys.: Condens. Matter 5, 2603 (1993). 

3 J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 

4 S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301, 
1348 (2003). 

5 J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, 
and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 (2004). 

8 Y. K. Kato et al, Science 306, 1910 (2004). 

7 J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, 
Phys. Rev. Lett. 94, 047204 (2005). 

8 N. A. Sinitsyn, E. M. Hankiewicz, W. Teizer, and J. Sinova, 
Phys. Rev. B 70, 081312(R) (2004). 

9 S.-Q. Shen, Phys. Rev. B 70, 081311(R) (2004). 

10 M.-C. Chang, Phys. Rev. B 71, 085315 (2005). 

11 A. V. Sh ytov, E. G M ishchenko, H.-A. Engel, and B. I. 
Halperin, cond-mat/0509702 

12 J. I. Inoue, G. E. W. Bauer, and L. W. Molenkamp, Phys. 
Rev. B 70, 041303 (2004). 

13 E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, Phys. 
Rev. Lett. 93, 226602 (2004). 

14 A. Khaetskii, cond-mat/0408136 

15 R. Raimondi and P. Schwab, Phys. Rev. B 71, 033311 
(2005). 

16 O. V. Dimitrova, Phys. Rev. B 71, 245327 (2005). 



17 O. Chalaev and D. Loss, Phys. Rev. B 71, 245318 (2005). 

18 K. Nomura, J. Sinova, N. A. Sinitsyn, and A. H. MacDon- 
ald, Phys. Rev. B 72, 165316 (2005). 

19 E. I. Rashba, Phys. Rev. B 70, 201309(R) (2004). 

20 S. Murakami, Phys. Rev. B 69, 241202(R) (2004). 

21 A. G. Mal'shukov, L. Y. Wang, C. S. Chu, and K. A. Chao, 
Phys. Rev. Lett. 95, 146601 (2005). 

22 A. Khaetskii, cond-mat/0510815 

23 B. K. Nikolic, S. Souma, L. P. Zarbo, and J. Sinova, Phys. 
Rev. Lett. 95, 046601 (2005). 

24 M. Onoda and N. Nagaosa, Phys. Rev. B 72, 081301(R) 
(2005). 

25 A. A. Burkov, A. S. Nunez, and A. H. MacDonald, Phys. 
Rev. B 70, 155308 (2004). 

26 J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

27 A. V. Zaitsev, Zh. Eksp. Teor. Fiz. 86, 1742 (1984) [Sov. 
Phys.-JETP 59, 1015 (1984)]. 

28 A. L. Shelankov, J. Low Temp. Phys. 60, 29 (1985). 

29 A. Millis, D. Rainer, and J. A. Sauls, Phys. Rev. B 38, 
4504 (1988). 

30 P. G. Silvestrov and E. G. Mishchenko, |cond-m at /0506516 

31 P. Schwab and R. Raimondi, Ann. Phys. (Leipzig) 12, 471 
(2003). 

32 P. Schwab and R. Raimondi, Eur. Phys. J. B 25, 483 
(2002). 



